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ABSTRACT 


The  paoer  presents  an  analysis  of  the 
stresses  and  strain*?  in  a fully  clastic,  rotating, 
annular  disc  that  has  initially  uniform  thickness 
and  is  made  of  a strain-hardening  material.  This 
analysis  is  based  on  Tresca1 s yield  condition  and 
the  associated  flow  rule,  and  assumes  that  the 
elastic  strains  may  be  neglected  in  comparison  to 
the  finite  plastic  strains  that  are  considered. 

The  bursting  speed  of  the  disc  is  expressed  in  the 
form  of  a definite  integral  which  involves  the 
strain-hardening  function  of  the  material.  In 
general,  this  integral  will  have  to  be  evaluated 
numerically,  but  analytical  evaluation  is  possible 
for  certain  strain-hardening  functions.  In  par- 
ticular, it  is  shown  that  for  linear  strain-harden- 
ing instability  can  occur  only  at  the  onset  of 
plastic  flow,  whereas  for  logarithmic  strain-hard- 
ening considerable  plastic  deformation  of  a stable 
character  may  occur  before  the  process  of  deforma- 
tion becomes  unstable  at  the  bursting  speed. 


♦The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N?onr-35301 * 

^Assistant  Professor  of  Applied  Mathematics. 
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NOTATIONS 

initial  inner  and  outer  radii  of  disc, 

instantaneous  (true)  rates  of  strain  in  the  radial, 
circumferential,  and  axial  directions, 

initial  uniform  thickness  of  disc, 

thickness  of  disc  during  plastic  flow, 

constants  in  logarithmic  stress-strain  relation 
Eq.  (16), 

initial  radial  distance  of  a particle, 

radial  distance  of  a particle  initially  at  rQ, 

nominal  stress  in  the  radial,  circumferential,  and 
axial  directions, 

critical  nominal  stress  (initial  yield  stress), 
radial  displacement  of  a particle, 
radial  velocity  of  a particle, 

nominal  strains  in  the  radial,  circumferential,  and 
axial  directions, 

critical  nominal  strain, 


density  of  disc  material, 

true  stress  in  the  radial,  circumferential,  and 
axial  directions, 

twice  the  critical  value  of  the  true  shearing  stress, 
angular  speed  of  the  rotating  disc. 
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(I)  Introduction 

The  analysis  of  stresses  and  strains  in  elastic  or 
partially  plastic  rotating  discs  has  been  discussed  repeatedly 
[1-6  ]*.  In  these  problems  tho  strains  are  of  the  order  of 
magnitude  of  elastic  strains  and  may  therefore  be  treated  as  in- 
finitesimal. As  a consequence  of  this,  the  equations  of  equili- 
brium may  be  satisfiel  for  the  undeformed  rather  than  tho  deformed 
disc.  Considerable  simplification  of  the  mathematical  wcrv  results 
from  this  approximation.  It  would  be  unrealistic,  however,  to 
treat  the  strains  as  infinitesimal  when  the  bursting  speed  of  a 
strain-hardening  plastic  disc  is  to  be  determined.  In  fact,  as  the 
angular  speed  of  such  a disc  is  gradually  increased,  considerable 
plastic  deformation  of  a stable  character  may  occur  before  the 
process  of  deforma  ion  finally  becomes  unstable  at  the  bursting 
speed. 

Since  the  prediction  of  the  bursting  speei  of  a strain- 
hardening  plastic  disc  requires  the  consideration  of  finite  plastic 
strains,  the  analyst  who  has  to  forego  the  simplification  resulting 
from  infinitesimal  strains  will  look  for  other  possibilities  of 
reducing  the  mathematical  work.  Thus,  Zaid  [7  ] uses  a deformation 
theory  of  plasticity  rather  than  a flow  theory,  in  spite  of  his 
realization  of  the  questionable  value  of  this  type  of  theory  [ 8 ] • 

The  alternative  analysis  offered  in  the  present  paper  has  recourse 
to  other  simplifying  assumptions.  Firstly,  it  is  assumed  that  the 
elastic  strains  are  negligible  in  comparison  with  the’  plastic  strains. 


♦Numbers  in  square  brackets  refer  to  the  bibliography  at  the  end  of 
the  pancr. 
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Secondly,  it  is  assumed  that-.  the  material  of  the  disc  obeys  Treses' s 
yield  condition  [9  land  the  associated  flow  rule  [ 10  ] . 

(II)  TRESCA'G  Yir.LD  CONDITION  AND  FLOW  RULE 

In  a rotationally  symmetric  stress  field  the  principal 

stresses  are  the  hoop  stress  0 , the  radial  stress  0 , and  the 

©’  r’ 

axial  stress  0 . The  principal  shearing  stresses  therefore  are 
z 

? ^ °0  “ °z  1 » 7 1 °r  " d?  I’  and  ^l°o  " °r  I • According  to  Tre  sea's 

yield  condition,  none  of  these  principal  shearing  stresses  can 
exceed  a critical  value  that  depends  on  the  stage  of  strain-harden- 
ing of  the  material.  Moreover,  for  plastic  flow  to  occur,  at  least 
one  of  the  principal  shearing  stresses  must  have  the  critical  value. 

If  only  one  of  the  principal  shearing  stresses  has  the 
critical  value,  the  flow  rule  associated  with  Trosca's  yield  condi- 
tion stipulates  that  the  Instantaneous  strain  rate  corresponds  to 
pure  shear  in  the  plane  of  maximum  shearing  stress.  The  sense  of 
this  shear  deformation  must  be  "appropriate",  i.c.,  it  must  corres- 
pond to  the  sense  of  the  maximum  shearing  stress. 

If  two  principal  shearing  stresses  attain  the  critical 
value,  the  flow  rule  admits  any  strain  rate  that  can  bo  considered 
as  resulting  from  the  superposition  of  appropriate  states  of  pure 
shear  in  the  two  planes  of  critical  shearing  stress.  It  should  be 
noted  that  in  both  cases  the  volume  is  preserved  during  plastic 
flow. 

To  express  this  yield  condition  and  flow  rule  analyt- 
ically, denote  by  0/°  the  critical  value-  of  the  shearing  stress  for 
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the  considered  state  of  strain-hardening,  and  by  Oq,  c , and  c„, 
the  instantaneous  rates  of  extension  in  the-  circumferential,  radial, 
and  axial  directions.  For  the  state  of  generalized  plane  stress 
occurring  in  a thin,  rotating  disc,  the  axial  strcs~  a is  zero. 
Moreover,  for  an  annular  disc  of  initially  uniform  thickness,  it 
may  tentatively  assumed  that  the  hoop  stress  d^  and  the  radial 
stress  cj  are  tensile  and  that  everywhere  exceeds  d . Under 
these  circumstances,  Trasca1 s yield  condition  and  tno  associated 
flow  rule  lead  to  the  folltwing  relations: 

°©  = 0 > °r  > 0,  dz  = °,  CD 

er  = 0,  Cq  = - cz  2 C.  (2) 

Since  dg  and  er  are  thus  knov/n  explicitly,  the  number  of  unknown 
functions  is  reduced  and  the  analysis  is  thereby  simplified. 

To  complete  the  description  of  the  assumed  mechanical 
behavior  of  the  disc  material.  It  is  necessary  to  formulate  a lav; 
of  strain-hardening.  For  simple  tension,  Let  the  stress-strain 
diagram  have  the  general  shape  of  the  line  OAB  in  Fig.  1.  The 
material  remains  rigid  until  the  tensile  stress  s reaches  the  value 

s*;  during  the  ensuing  plastic  flow,  the  stress  increases  with  the 

strain  e according  to 

s = fve).  (3) 


The  rate  of  strain-hardening  ds/de  has  the  initial  value  s*/e * 

(see  Fig.  1)  and  decreases  monetonically  with  increasing  strain. 

In  kq.  (j),  s and  e should  be  inluruiu. led  as  the  con- 


ventional stress  and  strain  computed 


• . -i  *. 

W X Oil 
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dimensions  cf  the  test  specimen.  On  the  other  hand,  the  quantities 
°Q’  °r»  °z’  d’  °©>  °r»  °z  introduced  above  snould  be  interpreted  as 
the  ^ruc  stresses  and  strain  rates;  they  are  defined  with  respect 
co  the  dimensions  in  the  considered  state  of  deformation. 

It  remains  to  discuss  the  application  of  the  strain- 
hardening law  (3)  to  the  states  of  stress  and  plastic  flow  occur- 
ring in  a rotating  disc.  The  state  of  plastic  flow  described  by 
(2)  is  independent  of  the  value  of  the  intermediate  principal  stress 
o • Moreover,  this  type  of  plastic  flow  is  possible  under  Trcsca* s 
flow  rule  when  the  state  of  stress  is  simple  tension  in  the  circum- 
ferential direction.  It  is  therefore  reasonable  to  assume  that 
Eq.  (3)  can  be  applied  to  the  problem  on  hand,  provided  that  s is 
interpreted  as  the  conventional  hoop  stress  s0  and  e as  the  conven- 
tional hoop  strain  e0. 

(Ill)  FULLY  PLASTIC  DISC 

Consider  a rotating  annular  disc  of  the  uniform  initial 
thickness  hQ,  the  initial  interior  radius  a0,  and  the  initial  ex- 
terior radius  b0.  From  the  elastic  stress  analysis  it  is  known 
that,  as  the  angular  speed  u of  the  disc  is  gradually  increased, 
the  yield  limit  is  first  reached  at  the  interior  surface.  At  some- 
what higher  speeds  there  will  be  an  inner  plastic  region  surrounded 
by  an  elastic  rc-gion.  Any  flow  that  may  occur  in  the  plastic  region 
would  have  to  satisfy  (2).  If  v = v(r)  is  the  distribution  of  the 
radial  velocity  at  a generic  instant  during  this  plastic  flow,  the 
radial  strain  rate  is  er  = 5v/5r.  Since  this  must  vanish  according 
to  (2),  the  radial  velocity  must  be  independent  of  r. 


Moreover, 
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since  all  elastic  strains  arc  neglected,  the  racial  di  spj acenent 
u at  the  elastic-plastic  interface  is  zero.  This  means  that  v and 
hence  the  circumferential  strain  rate  e^  = v/r  must  vanish  through- 
out the  plastic  region.  This  region  thus  remains  rigid  until  the 
elastic-plastic  interface  reaches  the  exterior  surface  of  the  disc. 
The  angular  speed  at  which  this  occurs  is  given  by  (see  Ref.  6, 
p.  103,  Eq.  9.53) 


(4) 


where  p is  the  density  of  the  disc  material.  The  fully  plastic 
stress  distri nation  at  the  onset  of  plastic  flow  is  readily  analyzed 
(see  Ref.  6,  p.  104,  Eq.  9.59  and  7ig.  9.5).  It  is  found  that,  at 
this  instant,  the  radial  stress  nowhere  reaches  the  yield  stress 
regardless  of  the  values  of  bQ  > a0  > 0.  Thus,  (1)  and  (2)  apply, 
initially  at  least,  during  the  ensuing  plastic  flow. 

Let  r0  denote  the  initial  radius  of  a particle  that  is 
found  at  the  radius  r when  the  speed  is  u>.  If  the  radial  displace- 
ment is  denoted  by  u,  then 


r - rQ  + u.  (5) 

Here,  r and  u may  be  considered  as  functions  of  the  independent 
variables  rc  and  u?.  (It  is  convenient  to  use  rather  than  to, 
since  the  sense  of  the  plastic  flow  is  Independent  of  the  sense  of 
rotation  or  the  sign  of  the  angular  speed.) 

As  has  been  shown  above,  the  first  Eq.  (2)  requires 
that  all  particles  have  the  same  radial  velocity  at  any  given  instant. 


-miL JStX  « ■— ■ ' — ■ J -■  ■■  - -- - 
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This  implies  that  the  radial  displacement  depends  only  on  w but- 
not  on  rQ.  Thus,  the  material  bounded  in  the  undefo^ned  state  by 
the  coaxial  cylinders  of  radii  rr  and  r0  + drQ  undergoes  the  radial 
displacement  u and  is  then  bounded  by  cylinders  of  radii  r = rQ  + u 
and  r + dr  = r + drQ  in  the  deformed  state.  Simultaneously,  the 
original  thickness  hQ  decreases  to  h.  Because  of  the  incompressi- 
bility of  the  material,  hQr0dr0  = hrdr  or,  since  dr  = dr0 

r.  r, 

h = lu  —2  = h. 


+ u 


The  true  radial  stress  dr  is  transmitted  across  a 


(6) 


section  that  is  nroportional  to  hr.  Since  hr  = hor0,  by  (6),  the 
conventional  and  true  radial  stresses  have  the  same  values  or  - sr. 
The  true  circumferential  stress  fig  is  transmitted  across  a section 


that  is  proportional  to  hdr;  the  corresponding  value  in  the  unde- 
formed state  is  h_dr  = h-dr, 

(J  c o 

stress  is  therefore  given  by 


formed  state  is  h^dr  = hrtdr.  The  conventional  circumferential 

(J  c o 


se  =do^  • 

“o 


(7) 


This  stress  is  relate  l to  the  conventional  circumferential  strain 


_ u 


(8) 


by  an  equation  of  the  form  (3). 

mhe  eouation  of  enuilibrium  in  the  deformed  state  is 


t’(rhd  ) ? 2 

= hog  - hpufr  . 


dr 


(9) 


With  reference  to  the  undeformed  state  this  equation  may  be  written 
as  follows: 
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5<roSr> 

arc 


se  - pu2r0(r0  4 U) 
f(JL)  - puTro(r0  + u). 


(10) 


Since  s vanishes  at  the  interior  and  exterior  surfaces  of  the  disc, 
i.e.,  for  r0  = a0  and  rQ  = b0,  the  integral  of  the  right-hand  side 
of  (10)  between  the  limits  aQ  and  bQ  must  vanish.  V.'ith  a = b0/aQ, 
r)  = u/aQ,  and  £ = r0/a0,  this  condition  yields 


na 

Ji  f(iA)<tt 

2 (a 3 - 1)  4 3n(a2  _ l)  ’ 


(ID 


For  a ’'■novn  strain-hardening  function  f and  given 

initial  dimensions,  the  right-hand  side  of  (11)  must  be  evaluate', 

analytically  or  numerically,  for  a set  of  values  of  r)  = u/a0.  Each 

such  evaluation  furnishes  one  point  of  the  plot  & vs.  u.  For 

q - 0,  in  particular,  l’(0)  = s*  and  the  integral  in  (11)  equals 

(a  - l)s*.  Equation  (11)  therefore  yields  the  correct  value  (4) 

for  the  angular  speed  at  which  plastic  flow  begins. 

o 

As  long  as  an  increase  in  w is  required  to  produce  an 

increase  in  u the  considered  plastic  flow  is  stable.  The  bursting 

2 

speed  <*>2  corresponds  to  the  maximum  of  the  curve  w vs.  u. 

The  preceding  analysis  is  based  on  the  assumption  that 

the  radial  stress  is  smaller  than  the  circumferential  stress  d-.. 

v y 

A formula  for  the  maximum  radial  strops  occurring  at  a given  speed 
u>  can  be  obtained  by  carrying  out  the  differ  ntiation  on  the  left- 
hand  side  of  (10)  and  equating  9sr/3rc  to  zero.  Thus, 
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p 

max  s = 5.  - nu  r fr  + nl.  f 1? ) 

where  rQ  is  now  the  initial  value  of  tho  radius  at  which  the  maximum 

of  <5,  occurs.  Equation  (12)  shows  that  oven  this  maximum  of  s is 

A r 

smaller  than  the  corresponding  value  of  s0.  How,  <3  = s and 

w r r 

0q  > s q,  by  Eq.  (7).  The  maximum  of  or  is  therefore  smaller  than 
the  corresponding  value  of  <Jq. 

(IV)  LINEAR  STRAIN -HARD^ - ING 
As  a rule,  the  integral  in  (11)  must  be  evaluated 
numerically,  but  strain-hardening  laws  can  be  devised  that  mabe 
analytical  evaluation  possible.  The  simplest  law  of  this  bind  Is 
represented  by  tho  line  OAC  of  Fig.  1.  With 


f(e)  = s*(l  + e/e*), 


(13) 


Eq.  (11)  yields 

2 2 

pa0u> 

6s* 


a - 1 4-  ( q/e  * ) log  a 
2(a3  - 1)  + 3q(a2  - 1) 


(1*+) 


Differentiation  of  (1*+)  with  respect  to  q shows  that  d(u^)/dq  = 0 
when 


e*  = 2Cq2  + q + i)  log  a 

3(a  - 1) 


05) 


Equation  (15)  establishes  tho  critical  rate  of  linear 
strain-hardening.  For  a given  value  of  a = b0/c  , the  strain-hr.rd- 
ening  parameter  e*  (see  Fig.  1)  must  not  exceed  the  value  (15)  if 
plastic  deformation  is  to  be  stable.  In  Fig.  2,  e*  as  given  by 
(15)  is  plotted  versus  a.  As  a — > 1,  i.e,,  as  the  width  of  the 
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annulus  tends  towards  zero,  the  critical  value  of  e*  tends  towards 
1.  This  result  is  familiar  from  Laszlo' s paper  fll  ]. 

The  sign  of  dCu^/dp  as  found  from  (l4)  does  not  change 
with  increasing  q . This  means  that  for  linear  strain-hardening  in- 
stability can  occur  only  at  the  onset  of  plastic  flow.  A disc  that 
starts  to  flow  in  a stable  manner  will  continue  to  do  so  as  the 
angular  speed  is  increased. 

(V)  LOGARITHMIC  STRAIN-HARDENING 
Another  strain-hardening  law  for  which  the  integral  in 
(11)  can  be  evaluated  analytically  is  given  by 

f (e  ) = k log  (t  + me),  (l6‘ 

where  k,  -t,  and  ra  are  constants.  These  may  bo  chosen  so  as  to  fit 
an  experimental  stress-strain  curve.  With  the  function  (16),  Eq . 


(11)  yields 


2.  2 


pa0tu 


(a-t  + mT))log(a^  + raq)  - (£  mq)log(-t  + mq)  - a-t  log  a 


ran 


2(a3  - 1)  + 3r)(a2  - 1) 


(17) 


Figure  3 shows  the  actual  stress-strain  curve*  for 
A-t  24S-T4  in  simple  tension  (full  line)  and  the  curve  obtained  from 
Eq.  (16)  with  k = 9,780  psi,  t = 77,  ra  = 8,060.  Figure  4 shows  the 
left-hand  side  of  (17)  vs.  q = u/a0  for  a = b0/ac  = 2.  The  bursting 
speed  w2  is  obtained  from  the  maximum  point  of  this  curve.  As  can 
be  seen,  the  disc  bursts  for  q - n2  = 0.2$.  Since  n = u/aQ,  and 


♦The  authors  arc  indebted  to  Professor  E.  D'Appolonia,  Carnegie 
Institute  of  Technology,  Pittsburgh,  Fa.,  for  the  use  of  this  curve, 
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e = u/rQ,  e = r)  aQ/rQ;  and  since  the  maximum  nominal  hoop  strain 
occurs  at  rQ  = a , the  disc  would  burst  at  a hoop  strain  of  251, 
From  Fig.  3*  it  can  be  seen  however,  that  this  hoop  strain  is  be- 
yond the  range  of  validity  of  the  experimental  stress-strain  curve. 
A more  realistic  picture  is  obtained  by  integrating  the  right-hand 
side  of  (11)  numerically,  using  the  experimental  stress-strain 
curve  in  ^ig.  3 (full  line)  ns  s = f(e).  Figure  5 shows  a plot  of 

p p 

pa  u /6s*  (where  s*  is  as  indicated  in  Fig.  3)  vs.  n,  fora  = 2. 

o 

Again,  the  bursting  speed,  is  obtained  from  the  maximum  point 

of  this  curve,  which  occurs  at  r]  = r\?  = 0.14,  and  thus  at  a maximum 
nominal  hoop  strain  of  ibt.  This  maximum  strain  occurs  before 
necking  takes  place  in  the  test  specimen,  and  thus  is  in  the  range 
of  validity  of  the  experimental  stress-strain  curve.  The  results 
just  obtained  illustrate  the  fact  that,  for  logaritlimic  strain- 
hardening  as  well  as  actual  strain-hardening  of  a similar  typo, 
there  is  a considerable  range  of  stable  plastic  deformation  before 
bursting  occurs. 
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